clc,clear,close all
load('pic1_2_cz.mat');
A=fspecial('gaussian');
i=imadjust(z,[0.2;1],[0.8;1],0.5);
imshow(i);
for n=1:2%两次高斯滤波
i=filter2(A,i);
end
for n=1:10
i=medfilt2(i);%十次中值滤波
end
imshow(i)
% figure;

a=edge(i,'prewitt');%使用prewitt算子进行边缘检测
figure;
% imshow(a);
